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ABSTRACT 

Context. Some photometric studies of extragalactic globular cluster (GC) systems using the optical and near-infrared 
colour combination have suggested the presence of a large fraction of intermediate-age (2-8 Gyrs) GCs. 
Aims. We investigate the age distributions of GC systems in 14 E/S0 galaxies. 

Methods. We carry out a differential comparison of the (g — z) vs. (g — K) two-colour diagrams for GC systems in the 
different galaxies in order to see whether there are indications of age differences. We also compare the different GC 
systems with a few simple stellar population models. 

Results. No significant difference is detected in the mean ages of GCs among elliptical galaxies. SO galaxies on the other 
hand, show evidence for younger GCs. Surprisingly, this appears to be driven by the more metal-poor clusters. The age 
distribution of GCs in NGC 4365 seems to be similar to that of other large ellipticals (e.g. NGC 4486, NGC 4649). Padova 
SSPs with recently released isochrones for old ages (14 Gyrs) show less of an offset with respect to the photometry than 
previously published models. 

Conclusions. We suggest that E type galaxies assembled most of their GCs in a shorter and earlier period than SO type 
galaxies. The latter galaxy type, seems to have a more extended period of GC formation/assembly. 

Key words, galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: star clusters 



1. Introduction 

Whenever a major star formation burst occurs, globular 
clusters (GCs ) are expected to form along in significant 
amou nts (e.g.|Ashman fc Zepf|1992 . Elmegreen fc Efremov 



1997j|Forbes et al.||1997b| |Brodie fc Strader||2006| . In the 
ACDM universe examples of such events that shape the 
evolution of galaxies are infall of gas into single dark mat- 
ter halos and (gas rich) mergers. The latter, are the main 
mechanism behind the hierarchi cal merging scenario (e.g. 
White|[l978l|Springel et alp005l ). 

GCs still remain the best approximations to simple stel- 
lar populations (SSP), despite the recent finding s of some 
clusters hosting multiple stellar populations (e.g. Piotto et 
al.|2007} in the Milky Way). It is difficult to resolve the in- 
dividual stars of a galaxy beyond the Local Group. On the 
other hand, GCs can be easily detected at distances as far as 
the Virgo and Fornax clusters. Hence, they have become a 
powerful tool in the study of galaxies in the local Universe. 
Observations of GCs translated into for example, ages and 
metallicities are now widely used to put constraints on th e 
formation and evolution galaxies ( Brodie fc Straderj|2006 1 . 



The importance of dating GCs stems from the possibil- 
ity of retracing back the major epochs of star formation in a 
galaxy. For extragalactic, unresolved GCs, only integrated 
properties can be measured. Optical colo urs suffer from 
the well known age/metallicity degeneracy dWortheyJ 19941 ). 
This degeneracy can in principle be lifted by combining op- 
tical and near- infr ared (NI R) colours (e. g. Puzia et al. 120021 
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|Hempel et al. 2003) [Larsen et al.||2005[ ). Generally studies 
using the optical/JNIK, colour combination have deri ved ages 



by mainly comparing GC colours with SSP models ( Hempel 
l et al.||2007a[ |Hempel et~aLl|2007b| |Kotulla et al.||2008[ 
|Brodie fc Strader|2006 and references therein) . The typical 
conclusion found from these studies is that these galaxies 
host mainly old clusters (with ages 10 Gyrs), although 
non-negligible fractions of younger clusters (—2 — 8 Gyrs) 
have been claimed in some cases: NGC 1316, NGC 4365 and 
NGC 5846. In NGC 1 316 the presence of inte rmediate-age 
(-2-3 Gyrs) clusters ( |Goudfrooij eTaL|2001b[ ) is consistent 
with its very young luminosity weighted age (Kuntschner 



2000 ). This implies that NGC 1316 has experienced a major 



merger event 3 Gyrs ago ( Goudfrooij et al.|2001a[ ). However, 
in NGC 4365 and NGC 5846 there is no evidence from inte- 



grated light (Yamada et al. 2006 and Thomas et al. 2005 
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respectively) that these galaxies have undergone recent star 
formation activity. 

A general characteristic of GC systems of massive 
galaxies is having their (optical) colours distributed in a 
bi(multi)-modal way. Spectroscopic studies have suggested 
that colour bimodalit y is fundamentall y due to metallic- 

as most GCs are 



to each other in the ACS images were flagged and if they 
appeared as one bigger source in the LIRIS frame they were 
also excluded. We stress that any results for NGC 4382 and 
NGC 4473 should be taken with caution as their K-band 
observations were obtained in highly extincted conditions. 



ity differences (e. g. |Strader et al. 
found to be old (Strader et al 



2007 



C( 



although age differences of ~2 Gyrs are still allowed within 
uncertainties. The colour distributions of GCs in NGC 4365 



ohen et al ||1998| , 3. Deriving ages of GCs 



and NGC 5846 d o not show the bimodality feature ( Brodie 
& Strader 2006[ ) as clearly, as other high mass galaxies. 
Instead, NGC 4365 has been reported to have a colour dis- 
tribution with a normal blue peak if compared to other 
galaxies but a broader red peak shifted toward s bluer (V-I) 



colours (Larsen et al. 



(2006) report bimodality for the (g 



2001 Larsen 



et al.|2005| ). |Peng et 

— z) colour distribut 



aution 

of this galaxy with a broad red peak. The colour distribu- 
tion of NGC 5846 is bimodal although the distribution of 
the central clusters shows a red pea k with a hint of som e 
additional intermediate colour GCs (|Forbes et al. 1997a). 



Adding the high fraction of intermediate-age GCs to this 
argument, suggests that the GC systems of these galaxies 
are not like other similar type ellipticals. 

The apparent inconsistent results of NGC 4365 and 
NGC 5846 with no sign of intermediate-age stellar popu- 
lations hosting a large fraction of intermediate-age GCs is 
intriguing. GC studies and integrated light studies of the 
stellar population content of a galaxy should yield compati- 
ble conclusions. With this in mind we have studied the ages 
of GC systems in 14 E/S0 galaxies through optical/NIR 
photometry. We compare differentially a large number of 
GC systems with homogeneous data. The GC systems are 
also compared with SSP models. 



2. Observations and data 

The homogeneous survey of extragalactic GCs that we car- 
ri ed out with optical and NIR photometry was presented 

In the same article the data 



Chies-Santos et al. (2010a 



reduction techniques were described in detail, here, we will 
only briefly mention them. We obtained K (K$, 2.2/im) 
band imaging with the LIRIS spectrographer and imager 
( Acosta Pulido et al.||2003| IManchado et al.||2004) at the 
William Herschel Telescope (Will) from 2007 to 2009. The 
K imaging was combined with archival F475W (~ Sloan g) 
and F850LP (~ Sloan z) images from the Advanced Camera 
for Surveys (ACS) on board the Hubble Space Telescope 
(HST). Our sample is a sub-sample of the S AURON sample 
( de Zeeuw et alj 2002 1 and comprises bright elliptical and 



lenticular galaxies (M B < 19) with (m - M) < 32. Each 
galaxy was imaged in the K-band for about 3.4 hours and 
reduced with LIRISDR and standard IRAF routines. The g 
and z images were reduced with the MULTIDRIZZLE task 
in the STSDAS package. 

The GCs were detected in the ACS images and trans- 
formed to the LIRIS coordinate system. Following this, 
aperture photometry and effective radii (R e s) were mea- 
sured for the cl uster candidate s, with PHOT (Stetson 1987) 
and ISHAPE (Larsen 1999), respectively. A final sam- 
ple of GCs was obtained applying the criteria: g < 23, 
0.5 < (g — z) < 2.0, 1 < i? off (pc) < 15. Finally, a careful vi- 
sual inspection was done, where obvious background galax- 
ies still left in the sample were excluded. Sources too close 



In this section we attempt to derive ages of GC systems 
employing two different methods. In Sect 3.1 we directly 
compare the GC systems with SSP models. In Sect 3.2 we 
look for age differences without making use of SSP models. 
The second method gives more promising results. 

3.1. Simple stellar population models 
3.1.1. Data-model comparison 

In Fig. [I] we plot 2-colour diagrams, (g — z) vs. (g — K) for 
the final sample of GCs, following the criteria outlined in 
Sect. 2. Padova SSP models are shown in this plot, with 
model sequences of constant age ranging from the lower 
right to the upper left of each panel (2, 3, 6, 14 Gyrs), see 
legend in the bottom panels. Lines of constant metallic- 
ity are also plotted. Metal poor clusters are expected to 
be located in the bottom left part of the diagram whereas 
metal rich ones should be located in the upper portion. 
Padova SSPs are retrieved from the CMP 2.2 input form 
(http://stev.oapd.inaf.it/cmd), with Marigo et al. (2008) 
isochrones, and are from now on referred to as Padova08 
SSPs. Compared with these models, the data appears to 
be consistent with the oldest track, 14 Gyrs. Nevertheless, 
the 6 Gyr intermediate-age track also seems to be a good 
match. Excluding NGC 4382 due to highly extincted data, 
qualitatively NGC 4406 seems to host GCs with perhaps 
younger ages than the rest of the sample, while NGC 4365 
does not appear to differ significantly from the other galax- 
ies in terms of the ages of its GCs. Differences and similar- 
ities of the GC systems are quantified in Sects. 3.1.2 and 
3.2. 

In Fig. 



2 Maraston (2005) and Chariot & Bruzual 



prep ) models are overpiotted in the 2-colour diagrams, as in 



Fig. [Tj The comparison of the data with these models (dif- 
ferent panels of Fig. [2]) would suggest younger ages than the 
comparison with the Padova08 models. Therefore, conclud- 
ing how old or young a GC is through direct comparison 
to SSP models does depend on the choice of the model. 
Nevertheless, ages estimates through Padova08 models are 
closer to spectroscopically derived ages. This is probably 
due to th e new AGB and T P-AGB treatment of |Marigo et] 



al. (2008[). Maraston (20051 uses a different prescription of 



the AGB and TP-AGB phases through the fuel consump- 
tion theorem. Chariot and Bruzual models use older ver- 
sions of the Padova tracks (e.g. Padova 1999) which include 
older calibrations of the AGB phase and do not include a 
detailed prescription of the TP-AGB phase. 



3.1.2. Blue and red peak ages 

In Fig. [I] the median values of the blue and red popula- 
tions illustrated with larger symbols are also shown. We 
have divided each GC system in a red (metal-rich) and a 
blue (metal-poor) sub-population. We have chosen to use 
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2.02.53.03.54.0 2.02.53.03.54.0 
(9-K) 

Fig. 1. 2-colour diagrams for the GC systems of the 14 
galaxies, (g — z) vs. (g — K). Padova08 SSP models are over- 
plotted with model sequences of constant age and metallic- 
ity according to the legend in the bottom panels. The ages 
increase from the lower right to the upper left. The metal- 
licity tracks increase from the lower left to the upper right. 
The blue and red larger symbols are the median of the blue 
and red sub-populations divided assuming a common di- 
vision value of (g — z) — 1.187. The green symbol in the 
panels corresponding to NGC 4365 and NGC 4278 indicate 
the median of the single populations as the (g — z) distri- 
bution of these galaxies was considered to be more likely 
unimodal. The red arrow on the bottom most right panel 
indicates a typical reddening vector. 



the median values as the observational scatter is colour de- 
pendant (with blue clusters having more spread) and the 
median is more robust against outliers. The value which in- 
dicates the division between blue a nd red clusters was d eter- 
mined by means of the KMM test ( | Ashman et al.|1994[ ). We 
fit two Gaussians to the (g — z) distributions of the galaxies. 
The value where the blue and red Gaussians intersect is the 
(g — z) division value for the sub-populations. We use the 
mean of this value for the two most cluster-rich galaxies, 
NGC 4486 and NGC 4649, (g - z) = 1.187. For NGC 4365 
and NGC 4278 the {g — z) distribution was found to be most 
likely unimodal. For this reason we also plot the median of 
the whole population with a different colour. The median of 
the blue sub-populations have a trend of falling on relative 
younger tracks than the median of the red sub-populations. 
For example the value which corresponds to the median of 
the red sub-population of NGC 4486 is located between the 
6 and 14Gyr SSP while the blue one falls between the 3 
and 6 Gyr SSP. 

In Fig. p5] a {g — K) vs (g — z) 2-colour diagram for 
the GCs of 12 galaxies combined, are shown excluding 
NGC 4382 and NGC 4473. The blue and red large sym- 
bols correspond to the median colour of the blue and red 
clusters respectively, assuming a fixed dividing colour at 
{g — z) = 1.187. This Fig. shows that the blue population is 
subject to larger scatter in the {g — K) direction than the 
red population. This is due to the fact that blue clusters 
are generally fainter in K, especially due to the selection in 
the g band. This is also seen in Fig. [T] The finding of on 
average younger blue GCs if compared to red ones in Figs. 
[T]and[3]is attributed to the direct comparison between data 
and SSP models. We stress that this should be interpreted 
with caution and verified with spectroscopy and/or with fu- 
ture SSPs. One might think that objects with large (g — k) 
colours and small error bars (see Fig. [IJ would be making 
the blue population artificially young. Part of these objects 
could be backgro und stars or foreground galaxies (see e.g. 
Puzia et al.||2004 ). We note that a careful visual inspection 



was perfor med removing obvious ba ckground galaxies and 
blends (see Chies-Santos et aT]|2010a ). There could still be 
some contaminants left, but only spectroscopic analysis of 
such objects would give more clear answers. 

In order to quantify the ages that best ma tch the GC 
systems we ran AnalySED ( Anders et al. |2004 ) on the me- 
dian magnitudes of the sub-populations, presented in Fig. 
[T] This code finds the best age and metallicity fits using 
model spectral energy distributions (SEDs). We used as in- 
put SEDs Padova08 SSPs integrated colours and kept the 
extinction fixed. For simplicity, Padova08 models starting 
at 1 Gyr were employed. At old ages, metallicity has a much 
stronger effect on the SEDs than does age, as seen in the 
2-colour diagrams shown before. Therefore, to get meaning- 
ful results it is important to have a grid in metallicity with 
sufficient resolution. For Padova08 SSPs one can download 
a grid of metallicities with a spacing that suits the user. 
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Fig. 2. 2-colour diagrams for the GC systems of the 14 galaxies, (g — K) vs. (g — z), as in Fig. |l Maraston (2005) with 
red horizontal branch at the two lowest metallicities (left panel) and Chariot & Bruzual (in prep I ( right panelj. 
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In Fig. [4] the best age and metallicity fits returned 
by AnalySED are shown. The error bars represent the 
standard error on the mean. We have used a metallicity 
grid with a spacing of A Z=0.0015. The ages derived by 
AnalySED vary between 1 and 14 Gyrs, i.e, the data is 
roughly consistent with the whole range of models. The 
upper limits in most cases however, appear to match our 



expectations (from spectroscopy, e.g. IStrader et al. 2005 



Cenarro et al. 2007 ) of the age of GCs in early type galax- 



ies consistent with ^10 Gyrs. We strongly caution against 
interpreting the ages derived in the present section as real 
ages, they only serve to give an impression of what one 
gets from a literal comparison of the data with SSP models. 
Mean ages are well below 10 Gyrs, formally the mean blue 
is ~3 Gyrs while the mean red is ^8 Gyrs. Note that this is 
the same qualitative conclusion reached with Figs, [l] and [3] 
the blue population is on average younger than the red one 
by using the Padova08 SSPs as a reference. The correspond- 
ing best metallicity fits arc shown in the right panel of Fig. 
[4} The median blue best fit is Z=0.00046 ([Fe/H] 1.61 




3.0 3.5 
(9-K) 

Fig. 3. The (g — K) vs (g — z) 2-colour diagram for the 
GCs of 12 galaxies combined, excluding NGC4382 and 
NGC 4473. The blue and red large symbols correspond to 
the median colour of the blue and red clusters respectively. 
The orange lines represent Padova08 SSPs of 2, 3, 6 and 14 
Gyrs, like in Fig. [I] 



dex) and the red is Z=0.00121([Fe/H] 1.19 dcxQ The 

metal-poor peak value is comparable, or somewhatlower 
than what is fou nd for typical GC sy stems of massive galax- 
ies. For example Cohen et al. ( |1998[) finds the blue and red 
metallicity peak values for the NGC 4486 GC system to be 
— 1.3 and — 0.7dex respectiv ely, with a dispersion of .3 dex. 
More recent measurements ( Brodie fc Straderj|2006 ) agree 
on the blue peak value, —1.3 but the red peak is signifi- 
cantly more metal-rich, —0.2 dex. 



3.2. Relative comparison without SSP models 

Until up to this point we have attempted to derive ages by 
direct data model comparison. It is well documented in the 
literature that a direct comparison of integrated GC colours 
with SSP models tends to yield ages that are sign ificantly 



younger than a Hubble time (e.g. Larsen et al.|2005[ ). While 
this effect is less pronounced in the most recent Padova SSP 
models compared to earlier models, we still find it difficult 
to believe the very young ages assigned to the blue GCs in 
particular (Fig. [4]). It is worth noting that the spacing be- 
tween isochrones of different ages remains small even when 
including near- infrared colours (Fig. [I]), so even small errors 
in the model colours will have a large effect on the ages. We 
expect that as the models continue to be refined, these off- 
sets will grow smaller. In the meantime, it may be fruitful to 
explore other ways of determining age differences between 
GC systems in different galaxies. We therefore now turn to 
a purely differential comparison, independent of any SSP 
models. 

The median of the age distribution of the GC system 
of NGC 4486 is reported to be 13 Gyrs with a dispersion of 



2 Gyrs, through spectroscopy ( |Cohen et al.|1998 ). The GCs 



of NGC 4649 are distributed in a similar way in the (<? — A") 
vs. (g — z) plane to those of NGC 4486. Being these two GC 



1 The approximation [Fe/H]~log(Z/Z©) and Z©=0.019 were 
used. 
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Fig. 4. The best age and metallicity fits returned by 
AnalySED with Padova08 SSPs for the median of the blue 
(metal poor) and red (metal rich) sub-populations indicated 
by blue and red symbols respectively. The median of a sin- 
gle population for NGC 4365 and NGC 4278 where KMM 
rejected the bimodal fit as opposed to the unimodal one is 
indicated by a green star. 



systems supposedly old, and the ones with more clusters 
(167 and 301 respectively) we take them as the fiducial old 
systems. We fit linear relations to the data in the (g — K) 
vs. (g — z) plane. These fits give more weight to the objects 
with smaller observational errors in both the abscissa and 
the ordinate. 

In Fig. [5] this best fit relation given by Eq. (1) is shown. 
In the same plot, the best fit for the blue and red sub- 
populations are also shown. These, are given by Eqs. (2) 
and (3) respectively. Eq. (2) is valid for (g — z) < 1.187 
and Eq. (3) for (g — z) > 1.187. One possible reason for 
the apparent discontinuity of the (g — z) vs. (g — K) rela- 
tions of the blue and red GCs is that the horizontal branch 
morphology is expected to change abruptly at intermediate 




2.0 



2.5 



3.0 

(g-K) 



3.5 



4.0 



Fig. 5. (g — K) vs. (g — z) for the joint GC systems of 
NGC 4486 and NGC 4649. The weighted best fit for the 
whole sample, given by Eq. (1) is over plotted in green. The 
best fit lines for the blue and red sub-populations, given by 
(2) and (3) are shown in blue and red respectively. 



colours-metallicities, as discussed in Sect. 4.2. A separation 
into blue and red best fit lines is taking this into account. 



(g -z) = 0.465 * {g - K) - 0.349 
(g - z) b = 0.260 *(g- K) b + 0.232 
(g - z) r = 0.340 * (g - K) r + 0.140 



(1) 

(2) 
(3) 



In Fig.[6]the best fit lines (1), (2) and (3) in the (g-K) 
vs. (g — z) plane for the joint GC systems of NGC 4486 
and NGC 4649 are over plotted for the different galaxies. 
Also plotted are the median colours for the blue and red 
sub-populations of the respective galaxies. This plot shows 
qualitatively how much and in what direction the mean of 
the sub-populations deviates from the lines in Fig. [5j For 
instance, the median of the blue and red populations of 
some GC systems, such as NGC 4406, deviate significantly 
to the right of the green line. 

In the following, we introduce a new method to derive 
relative ages of GC systems with the purpose of quantifying 
the magnitude and direction of the deviation of a GC from 
Eqs. (1) (2) and (3). In Sect. 3.1 we showed that relative to 
the SSP models, a younger cluster falls on the lower right of 
the (g — K) vs. (g — z) plane while an older cluster falls on 
the upper left of this diagram. By considering the green line 
as an SSP, with metallicities increasing from the lower left 
to the upper right, one can make statements about relative 
ages of GCs with respect to this line. We define 6, <5biue and 
<5 rc d as the projected distances of a GC from relations (1), 
(2) and (3) respectively. These, are given by: 



0.465 * (g -K)- [g - z) - 0.349 
\/0.465 2 + 1 

_ 0.260 * (g - K) - (g - z) + 0.230 
V0.260 2 + 1 



(4) 
(5) 
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Fig. 6. Same as Fig[l]but with the Padova08 SSPs replaced 
with the best fit relations given by Eqs. (1), (2) and (3) 
shown in Fig [5] over plotted. 



0.341 * (g - K) - (g - z) + 0.140 
V0.341 2 + 1 



(G) 



Fig. [7] shows a cartoon illustrating S for a cluster in the 
(g — K) vs. (g — z) plane. In this figure the location of 
younger and older clusters relative to this line as well as 
the location of the metal rich and metal poor clusters are 
plotted. The values of 5 are age tracers of GCs in this plane. 
An object falling to the right of the green line, will have a 
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Fig. 7. A cartoon illustrating the definition of S in the 
(g — K) vs. (g — z) diagram for one cluster (see text for 
details) and the location of the "younger" and "older" ob- 
jects relative to the green line. The location of metal-poor 
and metal-rich objects is also shown. 



positive 5 value. A positive 6 value means that this object 
has an age that is younger than the age that the green 
line traces (S = 0). In this sense a GC system whose mean 
difference is further away to the right of the S — value in 
the distributions are the systems with the best chances of 
hosting younger GCs. 

In Fig. [8] the histograms for the 6 values of the different 
galaxies are shown. By comparing the different S distri- 
butions one can check if there are age differences between 
the different distributions. It is readily seen that not all 
GC systems seem to be consistent with a 5=0. For exam- 
ple, note that the peak of the NGC 4406 S distribution is 

-0.2 whereas the peak of the NGC 4649 is 0.05. By 

taking the mean value of the distributions of S it is pos- 
sible to see how on average the different GC systems are 
younger (or older) from the joint GC system of NGC4486 
and NGC 4649. In Fig. [9] a dia gram shows the mean (and 
median) of 6 with the standard error on the mean for each 
GC system. The usual interpretation of any offset from 
(5=0 would be that there are (mean) age shifts. The larger 
shifts occur for 5 >0, i.e to younger ages. Therefore the 
GC system of NGC 4406 is on average younger than the 
GC system of NGC 4649 and NGC 4486. The best cases for 
hosting younger GCs are NGC 4382, NGC 4406, NGC 4526, 
NGC 4570 and NGC4660. All of them are SOs, except for 
NGC 4660 which is an E5, following classification from 
NED . However, thi s galax y is also quoted to be an E3/S0 
by Ferrarese et al. (2006) who report boxy isophotes and 
a faint blue structure of about 2.5' with two spiral arms 
within an ACS pointing. 

By performing the S analysis for the different GC sys- 
tems separate for blue and red clusters one can check which 
sub-population has a greater d shift. The means (and me- 
dians) of the <5 b i ue and c5 rec i distributions are plotted in Fig. 
[lOl The values of <5biuc and 6 rc d are based on the colour 
ranges where Eqs. (2) and (3) are valid. The blue popu- 
lations of NGC 4570, NGC 4526, NGC 4406 and NGC 4382 
have the greatest shift, <5bi U c 0.1. The red populations on 
the other hand, have an apparent smaller shift. Moreover, 
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Fig. 8. The 8 distribution: the difference between the GCs 
and the best fit line for the whole sample given by Eq. (1). 
The dashed line indicates the value (5 = 0. 



Fig. 9. The mean and median of the 8 distribution for the 
GCs of different galaxies, shown in Fig. 8} The error bars 
centred in the mean indicate the standard error on the 
mean of the quantity 8. The dotted line stands for (5 = 0. 
NGC 4382 and NGC 4473 are marked with yellow rectan- 
gles because their K-band data were obtained in highly ex- 
tinctcd conditions. 



the <5biuc offsets from 8\ 



blu 



in Fig. H] 



in Fig. 10 are very similar 



to the (5 offsets from <5 

It is important to check how much scatter around the 
vertical line of 8, Suuc and <5 ro d = might be attributed to 
zero point, photometric and aperture correction uncertain- 
ties. The accuracy of the LIRIS K-band photometric cali- 
bration is estimated to be ~ 0.03 — 0.05 mag, see Sect. 3.2 of 
Chies-Santos et al. (2010a). There are cases where the mean 
8 reaches values as large as three times the values of the 



uncertainty on the photometric calibration. Therefore only 
part of the scatter on Figs. [9| and [T0| could be attributed to 
uncertainties in the K-band photometric calibration. 

While measurements of the mean 8 parameter in dif- 
ferent galaxies reveal some interesting differences, we have 
not yet discussed how to interpret this in terms of actual 
age differences. To this e nd, we are guided by a comparison 
with SSP models. In Fig. 11 (g — K) vs. (g — z) 2-colour di- 
agrams are shown with lines of constant 8. Padova08 SSP 
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Fig. 10. Similar to Fig [9] but for the mean and median 
of the Sbiue and <5 re d distributions for the GCs of different 
galaxies. The error bars centred in the mean indicate the 
standard error on the mean of the quantities. The dotted 
line stands for the value $biuc and #red = 0. NGC 4382 and 
NGC 4473 are marked with yellow rectangles because their 
K-band data were obtained in highly extincted conditions. 



models are also overplotted, as in Fig. [T] Note that the 
best fit line (1), is roughly consistent with the 6Gyr track. 
Moreover the line of 5 ~ 0.1 is located between the 2 and 
3 Gyr old tracks. This means that an offset of 6 ~0.1 is 
roughly consistent with an upper limit of Aage~4 Gyrs. S 
values consistent with ^0.05 fall between the 3 and 6 Gyr 
old tracks, giving a maximum Aage~3 Gyrs. If we instead 
shift the model grid such that the 14 Gyr SSP model co- 
incides with the 6 = line, an offset of 6 = +0.1 would 
correspond to an age of about 6 Gyr, or a Aage of 8 Gyrs. 
It is thus clear that a translation of our 6 parameter to an 
actual age difference is quite uncertain. We stress that care 
must be taken on the interpretation of the <5-parameter as 
a pure proxy for age, due to the age-metallicity degeneracy. 




Fig. 11. (g — K) vs. (g — z) diagrams showing lines of con- 
stant S, according to the legend. Padova08 SSP models are 
over plotted as in Fig[l] The green line is the best fit line 
for the joint GC systems of NGC 4486 and NGC 4649, given 
by relation (1). 



For instance, the SSP tracks in the 2-colour diagram do 
not run parallel to the lines of constant Sbiue and 8 re d and 
the comparison between the two is metallicity dependant. 
Therefore it is difficult to even roughly estimate what ^biuc 
and (S re d mean in terms of Aage. 

One might ask whether differences in the 5 values could 
be related to metallicity rather than age effects due to the 
kno wn peak GC metallicity vs. galaxy luminosity relation 
(eg. Brodie & Strader 2006). This would happen because 
the fiducial best-fit lines Egs. (1), (2) and (3) are defined 
on two of the highest luminosity galaxies in our sample. 
We argue this does not appear to be the case. The rea- 
son is that the peak GC metallicity vs. galaxy luminosity 
relation does not translate into a peak GC metallicity - 
morphology relation. For example, the SO NGC 4406 which 
has a higher delta value has a Mb value comparable to that 
of galaxies with low delta values, eg. NGC 4365. Another 
example is NGC 4382 which has a high delta value and an 
M B value consistent with that of NGC 4649. Although the 
case for NGC 4382 is to be taken with caution due to the 
worse quality K-band data. Moreover, we have calculated 
the Sbiue for two clusters with the same age (10 Gyrs) and 
with host galaxy luminosities, Mb=-21(Z ~ —0.0 008) and 

M B =-19 GC (Z 0.0006), following Fig. 13 of |Peng et 

al 



(2006) 



This was done in order to test for the worst case 
scenario, the influence of the peak GC metallicity vs. galaxy 
luminosity relation on the measured 5 values. By making 
use of Padova08 SSPs we find that the difference between 
these <5biue values is 0.0064, i.e., much smaller than the ac- 
tual delta values. Therefore, the peak GC metallicity vs. 
galaxy luminosity relation has no significant impact on the 
d values over the luminosity range relevant for our sample. 



10 



Chies-Santos et al.: Age distributions of GCs in early-type galaxies 



Galaxy 


Morphological 


Morphological 




lype 


iype (plots J 


(9) 


<•?) 


(31 


Mriri 0077 
INOO 66 i I 


rjO/rj6 




Mpp /107Q 
IN 4Z / o 


rjl/ rjz 




NGC 4365 


E3 


E3 


NGC 4374 


El 


El 


ncc 438? 




S ARffl'l 


NPP AA0f\ 

1 N vi V. 'IIUU 


cm 


ou 


1 \ V I V ^±^± ( o 




pec 




rjU/rjl pec 




NGC 4526 


SABO (s) 


SAB(O) 


NGC 4552 


E0/E1 


E0/E1 


NGC 4570 


SO sp 


SO 


NGC 4621 


E5 


E5 


NGC 4649 


E2 


E2 


NGC 4660 


E 


E5 



Tabl e 1. G alaxy (1), morphological type from Emsellem 
et al. |2007 and NED (2) and morphological type used for 
Elgs.(l2jand[l3| 



3.3. GC ages correlated with host galaxy morphological type 



In Fig. 12 the mean (and median) of 5, Shine and r5 re d are 
plotted against the morphological type of the galaxy. The 
morphological classification of the galaxies is given in Table 
1. It is readily seen that the mean of S correlates with 
galaxy morphology. The S (and <5biue, ^rcd) distributions of 
the galaxies of later Hubble type (E5, SO's, SAB(O)'s) have 
on average larger positive offsets from S = (and <5biuc = 0, 
^rcd = 0) if compared to that of galaxies of earlier Hubble 
type (EO's). These offsets are interpreted as being due to on 
average younger GC ages. If the E5 type galaxy NGC 4660 
were to be classified as an SO, this correlation would be 
even stronger. To quantify the strength of this correlation 
we make use of the Spearman's rank correlation coefficient 
(p) which is a non-parametric measure of statistical depen- 
dence between two variables. If one of the variables is a 
perfect monotone of the other the quantity p will be equal 
to +1 for correlation and -1 for anticorrelation. Con sidering 
NGC 4660 as an E5(S0) in the relation of Fig. 



12 the val- 



ues of p are: 0.75(0.76) for the mean and 0.69(0.73) for 
the median. For 5 blue , p= 0.79(0.81), 0.81(0.83) and for 
^rcd, p=0. 35(0.41), 0.47(0.53) for the mean and median re- 
spectively. If one removes NGC 4382 and NGC 4473, due to 
highly extincted K-band photometry, the p values become 
0.67, 0.63. For r5 b i uc the p values are 0.74, 0.76 and for <5 rcd 
they are 0.41 and 0.65, for the mean and median respec- 
tively. The values of p for S (and <5biue) are close enough 
to 1 for the correlations between 5 (and <$biue) and galaxy 
morphology to be significant. 

From Fig. [IT] we have that an offset of S = 0.1 is roughly 
consistent with a Aage of ~ 4-8 Gyrs. If one consid ers the 
median age of ~13 Gyrs found by Cohen et al. ( 1998 1 for the 
GC system of NGC 4486, an E0, the average ages of SO GC 
systems are extrapolated to be ~5-9 Gyrs. We stress that 
this is a very rough estimate, given the zero-point problem 
of the models. 
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Fig. 12. The mean and median of 5, ^biuc and 5 rc d as a 
function of morphological type in the different panels. The 
solid line indicates the value where the mean of the GCs 
would not differ from the relations (1), (2) and (3). The 
arrows indicate the location of the values for NGC 4660 
when classified as SO instead of E5. 



4. Investigating effects other than age and 
metallicity in the 2-colour diagrams 



In the top panel of Fig. 13 the standard deviation (c) of the 
clusters around relation (1) along (g — K) in the (g — K) vs. 
(g — z) plane is shown as a function of morphological type. 
Because the photometric errors are dominated by the K- 
band photometry, the scatter in (g-K) will be very nearly 
proportional to the scatter in S if photometric errors are 
the only cause of offsets from the line defining the 5 param- 
eter. It is immediately seen that the scatter varies among 
different systems and that in most cases it is larger than 
the ler of ~ ±0.4 permitted ( Chies-Santos et al.| 2010a) 
from observational uncertainties. Analogously the scatter 
of <5bi uo and <5 rc d along {g — K) as function of galaxy mor- 
phology are shown in the middle and in the bottom panels 
respectively. It is readily seen that the scatter is greater for 
5 than for Jbiue and 8 rc< \ separately. For 6, the mean a is 
0.643 ± 0.007. The error quoted is the standard error on 
the mean. For f5 b i ue the mean is 0.518 ± 0.013 and for S rcd 
0.403 ±0.015. The scatter becomes smaller when looking at 
f5bi U e and r5 re d separately. In fact, it seems almost consistent 
with the observational errors. This suggests that a linear 
fit in the 2-colour (g — K) — (g — z) space is an over ap- 
proximation. By making the division in two linear regimes, 
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Fig. 13. The scatter represented by the standard deviation 
(a) for S, Shine and (5 roc i in the (g — K) direction as a function 
of morphological type. The dotted lines indicate the upper 
limit to the maximum scatter from observational errors. 




K<20 mag 



Fig. 14. The mean a (stars) an d the mean mad (squares) 
for the measurements in Fig. 13 for S (black symbols) , ^biuc 



(blue symbols) and 5 re d (red symbols). The error bars are 
the standard error on the mean. 



The a is consistent with ~0.4 for the first and drops below 
r~> 0.4 for the latter. However, a sample binned to brighter 
magnitudes will of course have a smaller mean observational 
error value. For K<19 mag this value is ~ 0.2. The a values 
for K<19 mag are much larger than this value for the single 
and for the red and blue sub-populations. Even though the 
mad values are closer to ~ 0.2 for the blue and red sub- 
populations, the observational errors are really a values and 
not mad values. Therefore there has to be something else 
responsible for this extra scatter. 

In this section we investigate which effects other than 
age and metallicity could be responsible for this extra scat- 
ter, and whether these are expected to contribute signifi- 
cantly to any scatter. 



it seems that almost all the scatter can be explained by 
observational errors. 

One might think that the relation shown in Fig 
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might be an artefact of our sky subtraction method through 
ellipse model fits (see Chies-Santos et al. 2010a). Even 
though the direction the morphology becomes more com- 
plicated is the same direction the ellipse model fits became 
more difficult, Fig. [T3| shows that there is no relation be- 
tween the scatter in S (Shiue and 5 re d) and the host galaxy 
morphological type. We therefore conclude that the ob- 
served trend of increasing 6 with Hubble type is real. 

In Fig. [14] a diagram with the mean scatter for all the 
GC sample, and for GCs brighter than K=20 and K=19 is 
shown. The scatter is now also represented by the median 
absolute deviation (mad, indicated by squares) in addition 
to the a (represented by stars). 

Both mean mad and a decrease when brighter mag- 
nitude intervals are considered. This is true both for the 
single population and when separating in the red and blue 
sub-populations. For the single population the values of a 
remain g reater than the observatio nal errors of ~ 0.4 (see 
Sect. 4 in Chies-Santos et al.|2010a[ ), even for a sample with 
K<19 mag. For the blue and red sub-populations, for the 
brighter magnitude samples, this is not the case any longer. 



4.1. a— enhanced vs. solar-scaled tracks 

What is the effect of a— enhancement on the int egrated 
colours of S SPs? For NGC 4472 and NGC 4486 GCs (Cohen 

q/Fe 



2003 and ICohen et al. 
d to be 



1998) ]q/Fe | estimates were 
+0.3. |Puzia"et al.|(2005|) finds an [a/Fe] 



et al 

derive^ 

~ +0.4 f or their sample of extr agalactic GCs. In the recent 
study of Sharina et al. (2010) it is found that the [a/Fe] 
ratios of GCs tend to be higher on average for giants if 
compared to dwarf galaxies. If a gE has recently accreted 
nearby dwarfs with its GCs, they should have therefore 
lower [q/ Fe] ratios than the GCs formed within the pro- 
togalaxy. Salaris & Cassisi (2007) discuss q-enhancement 
issues in the study of extragalactic GCs through broad- 
band photometry (see their Fig. 12) and conclude that the 
net effect of using an q-enhanced calibration in GCs with a 
solar metal mixture is to assign too young ages, especially 
at the reddest end. These age differences are of the order of 
~ 2 Gyrs. If on the other hand, one uses a solar scaled cali- 
bration to derive ages of q-enhanced GCs, the ages assigned 
will be too old. 

Recently the BaSTI group made available q— enhanced 
([q/Fe] ~ +0.35) isochrones and spectra. These were down- 
loaded from their website and used as input in GALEV to 
produce q-enhanced calibrated SSPs for the g, z and K 
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Fig. 15. A 2-colour diagram with a-enhanced and solar 
scaled SSPs for the ages of 2 and 10 Gyrs. The SSPs were 
produced in GALEV using spectra from the BaSTI group. 



colours. In Fig. 15 we plot 2 and lOGyr SSPs in a 2-colour 
diagram with solar-scaled and a-enhanced metal mixtures. 
The different tracks are so close together that it is very 
hard to tell them apart. The magnitude differences seen 
are too small to account for any noteworthy difference be- 
tween the two sets of models. The fact that Salaris fc Cassisil 
(2007) found a difference in comparing a-enhanced and 



solar-scaled SSPs using V, I and K magnitudes could be at- 
tributed to the filters themselves. Coelho et al. (20071 found 
that for the majority of colour combinations the effect of 
a is in the same direction of metallicity, but in some (e.g. 
(B-V)) the effect is opposite to metallicity. For (<? — z) and 
(g — K) the effect was unknown. We also tested the differ- 
ence between GALEV a— enhanced and solar scaled SSPs, 
produced using spectra from the BaSTI group with the V, 
I and K filters. The differences were slightly larger using 
the VIK filter combinations i nstea d of gzK, confirming the 
results of Salaris & Cassisi (2007). We conclude therefore 
that through GALEV a— enhanced SSPs using BaSTI spec- 
tra, [a/Fe] variations would not contribute significantly to 
the scatter (i.e., ~ 0.1 — 0.2 mag) in the 2-colour diagrams 
using the gzK filter combination, as done in this study. 

4.2. Horizontal-branch Morphology 



In Fig. 16 the [g — K) — (g — z) diagram for the joint 



NGC 4486 and NGC 4649 GC sample is shown with differ- 
ent models overplotted in different panels. Note the wavy 
feature the data presents around (g — K) ~ 3.2 and 
(g- z) ~ 1.2. A SpOT Teramo 1 4 Gyr SSP with a de 



tailed treatment of HB morphology ( Raimondo et aL|2005 ) 
is overplotted in the top panel. Even though this track does 
not fit well the data in the redder part of the diagram, it 
also presents the wavy feature seen in the data. Other SSP 
models that do not have a detailed prescription of the HB 
morphology do not show this behavior (middl e panel for 
Padova08 and bottom panel for |Maraston|2005 ) . As shown 
in Fig. [T3j the scatter in the S values is significantly reduced 
when going from a single linear relation between [g — z) and 
(g — K) to two linear relations. This separation into two lin- 




Fig. 16. (g - K) vs. (g - z) for GCs of NGC 4486 and 
NGC 4649, indicated by different symbols. Note the wavy 
feature the data presents around (g — K) ~ 3.2 and 
(g — z) ~ 1.2. The different panels have different SSP 
model over plotted: upper panel: SpOT- Teramo 14 Gyr SSP 
with a realistic treatm ent of horizontal branch morphology 
( Raimondo et al. 2005), middle panel: Padova08 and bottom 
panel: Maraston ( |2005 ). 



ear relations simulates the effect of varying HB morphology 
in the 2-colour diagrams. 

4.3. Stochastic effects 

What would fluctuations in the number of bright stars do to 
the integrated colours of extragalactic GCs? We have gen- 
erated synthetic integrated magnitudes and colours from 



models based on P adova isochrones (Marigo et al. 



Girardi et al. 2008) and a Monte Carlo method fol 



2008 



owing 



Santos & Frogel (1997) in order to assess the influence of 



stochastic effects on the integrated colours of the clusters. 
Similarly as the SSPs used in Sect. 3, Padova isochrones for 
g, z and K were retrieved from CMD 2.2. A set of isochrones 
combining 18 ages (2-14 Gyrs - log £=9.3 to 10.15 at 0.05 
steps) and 4 metallicities (Z=0.00019, 0.0006, 0.0019, 0.006) 
was employed. Integrated magnitudes and colors involving 



Chies-Santos et al.: Age distributions of GCs in early-type galaxies 



13 




Fig. 17. The synthetic (g — K) vs. (g — z) diagrams for 
10 6 stars, for four different mctallicities indicated by the 
different colours/symbols. Upper panel: for the ages of 2 
and 14 Gyrs. Bottom panel: for the age range of 2 to 14 Gyrs 
and with the GCs of NGC 4486 over plotted. 



g, z and K bands were computed by firstly distributing 
a prefixed number of stars on an isochrone of given age 
and metallicity according to a Salpeter initial mass func- 
tion (IMF) . The use of a Kroupa mass function would not 
make a significant difference. The stochastic nature of the 
IMF was simulated by randomly sampling stellar masses 
between 0.08 and 120 M weighted on the IMF until the 
preset number of stars was reached. One hundred SSPs were 
simulated this way with their integrated fluxes spread re- 
flecting the fluctuation in the number of bright stars. A 
GC with g=23 has ~ 6 x 10 5 M© at the distance of the 
Virgo cluster. This corresponds to ~ 10 6 stars, considering 
a Salpeter IMF. For simplicity, we have considered models 
with a fixed number of stars. We examine models of 10 6 
stars, since only the brightest GCs (g < 23) are in the 
sample. 

In Fig. 17 the (g — K)vs.{g — z) synthetic diagrams are 
Tr stars. The upper panel shows that the scatter 



shown for 

on the (g — K) direction due to stochastic effects produced 
by a 2 Gyr SSP is slightly larger than that corresponding 
to the older, 14 Gyr one. The stochastic fluctuations are 
slightly higher for the younger SSPs if compared to the 
older ones due to the fact that the main sequence turnoff 
point occurs at a higher mass for the 2 Gyr SSP. This gives 
fewer RGB/AGB stars for a given total mass due to the 
combination of shorter stellar lifetimes and the IMF. The 
bottom panel shows that even for the combined age range 



of 2-14 Gyr, there is not a good match between the SSPs 
and the NGC 4486 GCs. The stochastic effects seem to be 
negligible as a source of colour spreading for such popu- 
lous clusters. Most of the scatter due to stochastics is or- 
thogonal to the direction of varying age, not contributing 



significantly to the scatter as plotted in Figs. 13 and 14 



5. Discussion and implications for GC system 
assembly 

This study shows that the ages of the GC systems are cor- 
related with the morphology of the host galaxy. Galaxies 
of later type are found to have a GC system with an age 
shift towards younger relative values if compared to galaxies 
of earlier type. This finding can be crosschecked with host 
galaxy properties. Major star formation events should leave 
imprints both in the integrated light and in the GC system 
of a galaxy. Evidence of recent star formation in later-type 
galaxies i s also found in integrated light studies. Shapiro et 
al. (2010) combines SAURON integrated field spectroscopy 
of early-type galaxies with Spitzer/IRAC imaging to inves- 
tigate the presence of trace star formation in such systems. 
Connecting a detailed kinematic analysis given by the first 
data set with direct probes of star formation from Spitzer, 
these authors find that star formation in early-type galaxies 
happens exclusively in fast rotating systems. In Fig. [18] 5, 
<5biuc and <5 re( j are shown as function of Xr for the the differ- 
ent galax ies from Emsellem et al. (20071 and Cappellari et 



aL](2007). This para meter is a proxy for th e degree of ro 
tation of the systems ( Emsellem et al.|2007 ). Slow rotators 
have Xr below 0.1 whereas fast rotators above 0.1. Note 
that the GC systems consistent with younger ages, in our 
picture, fall mostly among the fast rotators, with A^ > 0.1, 
with a tendency for larger Xr values (> 0.45) for S and S^i uc . 
Slow rotators on the other hand seem to have a population 
of GCs consistent with 5 and S^iue < 0, implying old ages. 
This is suggestive of a more recent GC formation-connected 
to recent star formation in the host galaxy only occurring in 
fast rotating systems. The two outlying cases, with low Xr 
values but high 5 and (Jbiuc values are NGC 4406 (Xr < 0.1) 
and NGC 4382 (X R > 0.1). There is evidence that the fast 
rotator NGC 4382 could have an artificially l ow measured 
Xr in vi rtue of its significant inclination (see Emsellem et 
atlp007| . 

Star formation occurring later in time in a galaxy is 
expected to come from metal-enriched gas. Therefore, one 
would presume that the metal-rich clusters would be driv- 
ing this relation, yet this does not appear to be the case. 
The blue population appears to drive the relation between 
GC age and host galaxy morphological type. In this con- 
text, the blue GCs of the E5/S0s/SAB(0) are younger. This 
is surprising and contrary to the view th at the metal-poor 
GCs formed (exclusively) at high redshift ( |Brodie fc Strader 
2006). We have investigated whether this trend could be 



caused by systematics in the data such as seeing, aperture 
corrections, sky subtraction, etc., and conclude that this is 
not the case. There is, however a trend of anti-correlation 
with the number of clusters in a galaxy. However, this only 
reflects the fact that later type galaxies (SO's) have less 
clusters than earlier type ones (EO's). Furthermore, a trend 
between S values and the brightness of the cluster is not 
present. Here, the direction of S values shifts are interpreted 
as age shifts. However, given our limited understanding of 
(the necessary ingredients to construct) SSP models there 
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Fig. 18. S, <5biuc and <5 rc d vs. A#. parameter. The dashed 
vertical line corresponds to the division value between slow 
(Xji < 0.1) and fast rotators (A^ > 0.1). The solid line and 
the different symbols are as indicated in Fig. 12 



is still the possibility that this is not an age effect. It could 
be perhaps due to an unknown outcome of a certain group 
of stars in the clusters (e.g. HB, AGB). 

A possibility is that the blue population drives the re- 
lation between GC age and host galaxy morphological type 
as a result of minor mergers. There would exist two possi- 
ble ways for assembling these younger GCs. Dwarf galaxies 
have a more exte nded star formation history than galax- 
ies of larger mass ( |Tolstoy et al. 2009|) and contain almost 
exclusively metal-poor GCs (|l'orbes et al. 2000 1 . One pos- 
sibility is that these blue GCs had formed in the accreted 
dwarf galaxies ( Cote et al.|19 98 ) . Another possibility is that 
the GCs were formed during the merger event from the gas 
reservoir of accreted dwarfs, presumably with metal-poor 
gas. Evidence for this latter possibility is found in |Muratov| 
& Gnedin (2010). In this model for GC formation in the hi- 
erarchical framework of galaxy assembly, mergers of smaller 
hosts create exclusively blue clusters and mergers of more 
massive galaxies create both red and blue clusters. It re- 
mains uncertain, however, whether such effects would be 
expected to operate more prominently in SO-type galaxies 
compared to ellipticals. It is also unclear to what extent 
the age distributions of blue GCs in SO-type galaxies are 
directly related to the presence of recent star formation in 



the discs as revealed by SAURON. Unlike the blue GCs, 
the discs tend to be metal-rich, and spatially our blue GCs 
do not appear to be associated with the discs either. 

In the Milky Way, there is substantial evidence that 
part of its GC system is associated with accreted dwarf 
galaxies . The age-metallicity relation study of |Forbes &i| 
Bridges (20101 for galactic GCs suggests that there are 
two tracks of objects in the Galaxy. One with constant 
old ages ~ 13Gyrs, spanning a large metallicity range 
—2.2 <[Fe/H]< —0.6, possibly formed in situ. The sec- 
ond track, branches to younger ages and the objects appear 
clumped at a mean age value of ~ 10.5 Gyrs and intermedi- 
ate metallicity, [Fe/H] ~-1.3. This latter track seems to be 
dominated by GCs associated with dwarf galaxies accreted 
by the Galaxy. A 5 ~0. 05-0.1 corresponds in very rough 
terms to an age shift of ~2-8Gyrs. Such an age shift is 
therefore comparable to the difference in age between GCs 
formed in situ and GCs accreted from dwarf galaxies in 
the Milky Way. Moreover the [Fe/H] value (—1.3) of this 
clump of objects, associated to dwarf galaxies, is closer to 
the metal-poor peak metallicity value than to the metal- 
rich one of GC systems of early-type galaxies. For exam- 
ple, the metal-poor peak of NGC 4406 is at [Fe/ H]= -1.23 
whereas the metal-rich one at [Fe/H]= —0.56 (Brodie & 
Stradcr 2006, based on the empirical relation of |Peng ~ 



al.||2006| ). Therefore if this younger clump of objects of the 
Galaxy was to be placed in NGC 4406, for example, it would 
belong to its metal-poor population. This is suggestive of 
younger GCs possibly belonging to the metal-poor popu- 
lation of a galaxy. In this context the age and metallicity 
of the blue population of clusters are related to the objects 
associated to dwarf galaxies. 

There is also a trend of later type galaxies exhibiting 
a greater po sitiv e <5 ro( j shift than earlier type ones (bottom 
panel of Fig. 12 ), although with much more scatter than for 
^biue- The maximum shift occurs for E5 galaxies. A great 
5 re d shift could be indicative of the host galaxy having ex- 
perienced a major merger. From the Ashman & Zepf (Tl992) 
scenario, formation of metal-rich GCs would occur during 
the merger of disc galaxies. Alternatively, it is tempting 
to link the large S YC d values to the amount o f "disc" for- 



mation tha t happened in a particular galaxy (Kuntschner 
et al.||2006|). For example, the fast rotators with large d re d 



values, NGC 4660, NGC 4473 and NGC 4526 have a promi- 
nent Mgfr enhanced, metal-enriched central disc-like struc- 
ture ( |Kuntschner et al.||2006||Kuntschner et al.|[2010| . 

For all this to get beyond speculation, further investiga- 
tion is required. Increasing the number of measured <5's for 
different GC system would perhaps strengthen (or weaken) 
the relation between GC system age and morphological 
galaxy type. Spectroscopy of GCs in later type early-type 
galaxies (SO's, perhaps E5's) might give it further confirma- 
tion. Moreover, cosmological simulations up to GC system 
formation and assembly could give further constraints on 
the nature of the mechanisms that shape the age distribu- 
tions of GCs in present day early-type galaxies. 

This study also shows that the age distribution of the 
GC system of NGC 4365 is not different from other similar 
type galaxies (e.g. NGC 4486 and NGC 4649). The optical 
colour distribution of this GC system is indeed different 
from others. This does not imply however, significant dif- 
ferences in the ages, but rather in the metallicity distribu- 
tions. 
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6. Summary and Conclusions 

Using multi-band g, z and the newly derived K-band pho- 
tometry we have compared differentially the GC systems 
in a sample of 14 E/SO galaxies. We have also studied the 
effect of different parameters in various SSP models. The 
results are outlined below. 

1. The most recent Padova SSPs for old ages show less 
of an offset with respect to the photometry than previ- 
ously published models. However, a formal fit still yields 
intermediate ages. 

2. The age distribution of the GC system of NGC 4365 ap- 
pears similar to that of other large ellipticals in the sam- 
ple. Although its colour distribution is different from 
other systems in (g — z), this difference is orthogonal to 
the direction of decreasing/increasing age shift in the 
(.9 - K )~(g - z ) plane. 

3. When the Padova08 SSPs integrated colours are used as 
SEDs the fits of the median ages of the blue and red sub- 
populations with AnalySED are found to be consistent 
with a single age, within one sigma. This makes them 
consistent with a great range of age values (from 1 to 
14 Gyrs), with the exception of some SO galaxies. 

4. A differential comparison without using SSP models is 
performed. This comparison is done in the direction of 
age shift, traced by the quantity 5. A correlation be- 
tween galaxy morphological type and the mean delta is 
found. Galaxies of earlier type host genuinely old GCs 
and galaxies of later type, younger on average clusters. 
The blue population appears to be driving this relation. 
The relation for the red population is not as clear as for 
the blue one. 

5. The (g — K) vs. (z — K) parameter space is analysed and 
it is concluded that there is a non-linear feature around 
(g — K) <~ 3.2 and (z — K) ~ 2. This feature is suggested 
to be due to the HB morphology. It is also present in 
SpOT-Teramo SSPs which have a detailed modelling of 
this stellar evolution phase. 

6. The influence of stochastic effects is investigated to ad- 
dress the spread seen in the 2-colour planes. It is con- 
cluded that the observational uncertainties are larger 
than the influence of stochastic effects for the cluster 
mass relevant to our study. In fact, the observational 
uncertainties can account for most of the observed scat- 
ter when considering blue and red clusters as separate 
populations. 

7. We find that there is residual scatter in the delta param- 
eter at bright magnitudes that cannot be accounted for 
by observational errors. This may indicate an intrinsic 
scatter in the GC ages of up to several Gyr. 
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